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Abstract 


The  effects  of  transmission  path  (and  environment)  on  the  characteristics  of  a 
transient  signal  due  to  HPM  sources  as  it  passes  through  electronic  systems  are  of 
importance  in  discovering  how  the  digital  circuits  might  be  altered,  since  such  signals  are 
unexpected  and  contain  spurious  electromagnetic  energy.  Theses  effects  can  be  studied 
by  numerically  modeling  the  propagation  of  the  transient  signals  through  electronic 
systems  such  as  a  mock  personal  computer.  For  estimating  the  wideband  system  response 
most  efficiently,  a  preliminary  analysis  is  presented  on  the  use  of  extrapolation 
techniques  to  extrapolate  in  both  time  and  frequency  domains  simultaneously  using 
Orthonormal  associate  hermite  polynomials. 
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Chapter  1 
Introduction 


The  influence  of  the  transmission  path  and  environment  on  the  characteristics  of  a 
transient  signal  as  it  passes  through  an  electronic  system  is  of  importance  in  discovering 
how  the  operation  of  digital  circuits  might  be  altered  as  a  result  of  unexpected  or  spurious 
electromagnetic  energy  being  incident  upon,  and  propagating  through,  an  electronic 
system.  If  a  signal  of  a  specified  form  enters  a  system  at  a  given  point  it  is  useful  to  know 
the  salient  features  of  the  resulting  signal  that  reaches  some  location  within  the  system, 
where  a  susceptible  digital  circuit  could  be  located.  Moreover,  it  is  useful  to  know  if  any 
of  the  features  of  the  induced  signal  are  primarily  those  peculiar  to  the  entering  signal,  or 
if  they  are  more  influenced  by  the  properties  of  the  transmission  path  and/or  environment. 
Addressing  these  issues  should  help  one  gain  an  appreciation  for  the  nature  of  a  spurious 
signal  arriving  at  the  input  of  a  digital  circuit  embedded  deep  within  a  complex  system 
and  perhaps  this  information  can  be  used  for  designing  and  building  better  high 
performance  systems. 

Firstly,  the  influence  of  the  transmission  path  (and  environment)  on  the  shape  and 
magnitude  of  a  transient  signal  is  considered  by  numerically  modeling  the  signal 
propagation  through  a  variety  of  structures  that  contain  certain  features  present  in  the 
real-world  systems.  Structures  that  electromagnetically  replicate  only  the  important 
details  of  the  component  configurations  found  in  typical  real-world  systems  are 
considered.  Examples  of  such  systems  include  aircrafts,  missiles,  ground  vehicles, 
personal  computers,  and  electronic  test  equipment. 


With  the  powerful  Maxwell  equations  solvers  available  today  to  perform 
eleetromagnetic  simulation  it  is  reasonable  to  assume  they  can  be  useful  in 
characterization  of  such  effects.  However,  given  the  complexity  of  real-world  systems 
this  still  remains  a  daunting  task,  especially  when  wideband  characterization  is  needed. 
For  example,  the  use  of  the  popular  FDTD  method  for  modeling  electronics  inside 
enclosures  imposes  serious  memory,  time,  and  accuracy  constraints  in  generating  the 
late-time  response.  Also,  the  use  of  MOM  or  FEM  when  generating  the  high-frequency 
response  can  be  computationally  overwhelming.  Recently,  researchers  have  reported  on 
wideband  extrapolation  techniques  which  can  potentially  overcome  some  of  these 
limitations  by  using  only  early-time  and  low-frequency  data.  By  using  only  the  early  time 
and  low  frequency  data,  which  contain  mutually  complementary  information,  one  can 
simultaneously  extrapolate  the  entire  system  response  in  both  domains. 

Finally,  we  shall  examine  how  well  these  extrapolation  techniques  work  for 
computing  the  signal  that  penetrates  an  enclosure  due  to  an  exterior  transient  source  such 
as  a  HPM  source.  The  FDTD  method  will  be  used  to  obtain  the  early-time  response  and  a 
triangular  surface  patch  code  (or  some  other  MOM-based  code)  will  be  used  to  obtain  the 
low-frequency  response.  Hence,  the  numerical  models  will  represent  different 
discretization  of  the  enclosures.  To  further  reduce  complexity,  only  the  features  deemed 
electromagnetically  essential  and  tractable  will  be  included  in  the  models.  Attempts  will 
be  made  to  ascertain  the  accuracy  of  the  approach  as  well  as  the  CPU  time  savings 
measured  against  using  either  the  FDTD  or  MOM  tools  separately. 
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Chapter  2 

Transient  Signal  Penetration  Into  Electronic 


Systems 

I.  Introduction 

Electromagnetic  interference  has  become  a  limiting  factor  in  digital  computer 
performance.  The  effects  of  electromagnetic  waves  on  highly  sensitive  digital  circuits 
composed  of  minute  circuit  elements  are  very  complex  and  difficult  to  understand. 
Simple  mismatches,  slots,  or  discontinuities  can  cause  coupling  in  the  printed  circuit 
boards.  Metal  traces  both  with  and  without  bends  are  likely  to  have  reactive  impedances 
that  degrade  the  performance  of  complex,  high-speed  circuitry  surrounding  them. 
Knowledge  of  such  electromagnetic  wave  effects  can  be  very  helpful  for  designing  and 
building  better  high  performance  systems.  Such  knowledge  can  also  be  helpful  in  saving 
human  lives  as  the  need  for  bloodshed  is  reduced  if  enemy  electronic  systems  can  be 
adversely  affected  through  proper  exposure  to  electromagnetic  interference.  Likewise, 
human  lives  can  be  saved  by  designing  systems  harden  against  such  threats. 

The  goal  of  this  chapter  is  to  predict  the  properties  of  a  transient  signal  as  it  makes 
its  way  through  a  propagation  path  from  the  exterior  of  a  structure  or  substrate  to  the 
terminals  of  a  deeply  embedded  digital  system  or  circuit  and  to  catalog  the  properties  of 
signals  that  might  be  expected  to  reach  deep  inside  digital  systems.  Firstly,  a  widely  used 
electronic  system  —  a  personal  computer  (PC)  is  studied  by  electromagnetically 
replicating  only  the  most  important  details  of  its  component  configuration  in  FDTD  using 


specialized  procedures  to  model  UPML,  thin  slots,  and  wire  antennas  as  described  in 
in  the  literature.  The  effects  of  the  transient  signal  propagation  through  the  modeled 
electronic  system  using  a  thin-wire  antenna  with  a  transmission  line  feed  to  excite  the  PC 
are  then  viewed  for  various  time-varying  waveform  excitations.  The  concentration  of  the 
transient  fields  at  various  locations  in  the  model  is  studied  to  access  more  information 
about  the  electromagnetic  wave  effects  at  those  locations.  Moving  a  step  forward  to 
further  analyze  the  effects  of  the  propagation  path  on  the  penetration  of  signals  into  the 
electronic  system,  a  bent  wire  structure  is  modeled  to  electromagnetically  replicate  a 
power  supply  cord  coming  from  the  outside  the  system  and  terminating  at  the  power 
supply  unit  and  the  disk  drives.  The  bent  wire  structures  terminating  at  the  power  supply 
unit  and  the  CD  drive  are  modeled  very  close  to  one  another  to  study  the  effects  of  the 
propagation  path  and  the  coupling  effects  of  closely  spaced  wire  bundles  in  real  PCs. 

11.  Computational  Model 

Consider  a  widely  used  electronic  system  like  the  small,  shoe-box  sized  Shuttle 
PC  shown  in  Fig.  2-1.  The  Shuttle  PC  was  considered  because  of  its  features,  which 
resemble  those  of  a  typical  desktop  PC,  and  its  compact  size.  The  features  of  this  PC  are 
closely  studied  in  order  to  electromagnetically  replicate  only  the  important  details  of  the 
component  configuration  found  in  typical  real-world  PCs.  The  electromagnetic  replica 
called  the  mock  PC  consists  of  conducting  walls  with  slots  for  air  vents  on  both  sides  and 
on  the  back,  hollow  PEC  cavities  for  the  various  devices  like  floppy,  CD,  and  hard  drive, 
and  the  power  supply  unit  as  shown  in  Fig.  2-2. 
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(a)  The  front  view 


(b)  The  rear  view 


(e)  The  side  view  showing  its  internal  struetural  eontent. 


Fig.  2-1.  Pietures  of  (a)  the  front  view,  (b)  the  rear  view,  and  (e)  the  side  view  showing  its 

internal  struetural  eontent,  of  the  Shuttle  PC 


The  eomputational  domain  for  the  FDTD  model  of  the  moek  PC  is  first 


terminated  with  PEC  walls  on  all  six  sides,  making  it  a  reetangular  PEC  box.  To  simulate 


an  open  spaee  region,  UPME  of  thiekness  15  eells  is  laid  on  the  surrounding  five  walls. 
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The  rectangular  chassis  of  the  mock  PC,  which  is  made  of  PEC,  is  then  modeled  on 

an  infinite  ground  plane.  The  dimensions  of  the  mock  PC  model  shown  in  Fig.  2-4  are  as 

mentioned  below: 

•  The  chassis  of  the  mock  PC  is  a  rectangular  PEC  box  of  dimensions 
30cm  X  20cm  X  20cm,  and  its  mesh  dimensions  are  60Ax40Ax40A  where 
Ax  =  Aj  =  Az  =  A  =  0.005  mand  At  =  0.0083  ns  . 

•  It  consists  of  the  following: 

•  CD-ROM  drive  of  dimension  32  A  x  30A  x  lOA 

•  Floppy  drive  of  dimension  32  A  x  20A  x  lOA 

•  Hard  disk  of  dimension  26 A  x  5A  x  30A 

•  Heat  sink  of  dimension  15Axl5AxlOA  located  lOA  above  the  ground 
plane 

•  The  power  supply  unit  of  dimension  1 5 A  x  1 OA  x  1 8 A 
These  are  modeled  as  crude,  hollow  blocks  of  PEC 

•  Illuminated  by  a  thin- wire  monopole  antenna  of  length  /  =  20  A  and  radius  a  = 

0.003m  located  15A  from  the  back  of  the  mock  PC  along  its  center  line  -see  Fig.  2-2. 

•  The  wire  antenna  is  fed  with  a  transmission  line  of  characteristic  impedance 
Zq  =  50Q  .  An  incident  waveform  is  introduced  in  the  transmission  line  using  a  one¬ 
way  injector  as  described  in  Chapter  5. 

•  There  are  vent  slots  of  width  0.5 A  on  both  sides  of  the  mock  PC  to  replicate  the  air 
vents  in  the  Shuttle  PC.  These  thin  slots  are  modeled  using  the  ETSF  mentioned  in 
Chapter  3. 
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•  As  indicated  in  Fig.  2-2  (c),  the  vent  slots  at  the  rear  of  the  moek  PC  for  the  heat 
sink  fan  are  of  thiekness  A  and  all  the  other  slots  on  the  power  supply  unit  are  of 
width 0.5A .  Here  again  the  sub-cellular  slots  are  modeled  using  the  ETSF  described 
in  Chapter  4. 
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(a)  The  over  all  eomputational  model 


Front  View 


Solid  Metal  Face  not  Shown 


Ax = Ay  =  Az  =  0.005m 

(b)  The  front  view 


Rear  View 


Ax=Ay  =  Az=0.005m 


(b)  The  rear  view 


Fig.  2-2.  (a)  The  over  all  computational  model,  (b)  the  front  view,  and  (c)  the  rear  view  of 

the  mock  PC 
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(a)  The  differentiated  Gaussian  pulse 


(b)  The  frequency  spectrum  of  the  differentiated  Gaussian  pulse 


Fig.  2-3.  (a)  The  differentiated  Gaussian  pulse  and  (b)  the  frequency  spectrum  of  the 
differentiated  Gaussian  pulse  introduced  into  the  transmission  line  feed  of  the  thin-wire 
antenna,  used  to  illuminate  the  mock  PC  from  behind 


III.  Transient  Signal  Penetration 
In  the  analysis  of  transient  signal  penetration  into  the  mock  PC,  the 
aforementioned  thin-wire  monopole  antenna  is  excited  at  is  base  by  the  differentiated 
Gaussian  pulse  of  Fig.  2-3.  The  cross-sectional  plane  bisecting  the  y-dimension  of  the 
mock  PC  in  Fig.  2-2  is  shown  in  Fig.  2-4  so  that  the  monopole  antenna  may  be  seen. 
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Fig.  2-4.  The  cross-sectional  plane  bisecting  the  y-dimension  of  the  mock  PC 


Viewing  the  fields  in  the  plane  bisecting  the  y-dimension  (see  Fig.  2-4),  one 
should  observe  zero  fields  inside  the  heat  sink,  the  CD-ROM  drive,  and  the  floppy  drive 
at  all  times,  as  shown  in  Fig.  2-5. 
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i  coordinate 


(a)  Time  step  =110 


(b)  Time  step  =  600 

Fig.  2-5.  The  field  for  transient  signal  penetration  through  narrow  slots  in  the  eross- 
sectional  plane  biseeting  the  y-dimension  for  (a)  time  step  =110  and  (b)  time  step  =  600 


From  Fig.  2-5  (a)  at  time  step  =110,  one  ean  see  that  there  is  a  significant  amount  of 
signal  penetration  into  the  mock  PC  through  the  vent  slots  in  the  back  panel  due  to  the 
near  proximity  of  the  monopole.  The  signals  thus  penetrating  the  system  get  trapped  in 
the  rectangular  chassis  and  bounce  back  and  forth  between  the  PEC  walls  of  the  PC  and 
the  crude  PEC  models  of  the  heat  sink,  CD-ROM,  floppy  drive,  and  hard  drive.  Eig.  2-5 
(b),  at  time  step  600,  shows  that  the  transient  signals  are  still  significant  in  the  PC  long 
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after  the  souree  has  died.  The  existenee  of  strong  fields  below  the  heat  sink  and  near 
the  sharp  corners  of  the  floppy  drive  and  the  heat  sink  might  be  a  likely  cause  of 
electromagnetic  disturbances  that  interrupt,  obstruct,  or  otherwise  degrade  or  limit  the 
effective  performance  of  electronics  in  the  PC. 


Fig.  2-6.  The  cross-sectional  plane  bisecting  the  z-dimension  of  the  mock  PC 
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(a)  Time  step  =120 


(b)  Time  step  =  600 


Fig.  2-7.  The  field  for  transient  signal  penetration  through  narrow  slots  in  the  eross- 
seetional  plane  biseeting  the  z-dimension  of  the  moek  PC  for  (a)  time  step  =  120,  (b)  time 

step  =  600 


Again,  in  the  eross-seetional  plane  biseeting  the  z-dimension  of  the  moek  PC  for 
Fig.  2-7  there  are  zero  fields  at  the  base  of  the  floppy  drive  at  all  times.  When  there  are 
strong  fields  near  the  monopole  at  time  step  =  120,  one  ean  see  a  signifieant  amount  of 
signal  penetration  into  the  moek  PC  through  the  vent  slots  in  its  panel  nearest  the 
monopole.  When  there  are  no  fields  near  the  monopole  at  time  step  =  600,  one  ean  see 


24 


the  fields  lingering  inside  the  moek  PC  long  after  the  source  has  died  out.  There  is 
very  little  field  penetration  into  the  power  supply  unit  at  time  step  =  120  and  these  fields 
become  significant  as  time  goes  on. 


Fig.  2-8.  The  field  Ey ,  which  enters  the  mock  PC  from  the  rear  panel,  showing  the 
transient  signal  penetration  through  the  narrow  slots  at  time  step  =  180. 


To  further  demonstrate  the  significant  transient  signal  penetration  into  the  system 
through  the  vent  slots  one  can  view  the  fields  penetrating  the  system  from  the  slots  at  the 
rear  end  of  the  PC  as  shown  in  Fig.  2-8.  Similarly,  the  very  minor  transient  signal 
dissipation  from  the  air  vent  slots  located  on  either  side  of  the  mock  PC  can  be  viewed  as 
shown  in  Fig.  2-9  (b).  The  field  ,  on  the  right  hand  side  panel  of  the  PC  at  time  step  = 
600,  also  illustrates  the  fields  lingering  in  the  PC  long  after  the  source  has  died  out. 
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Ax = Ay = Az  =  0.005m 

(a) 


I 
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(b)  Time  step  =  600 


Fig.  2-9.  (a)  The  right  hand  side  panel  of  the  mock  PC  shown  in  Fig.  2-4.  (b)  The  field 
on  the  right  hand  side  panel  of  the  mock  PC  shown  in  (a) 


To  obtain  a  better  understanding  of  the  electromagnetic  waves  entering  an 
electronic  system  from  an  outside  source,  the  field  intensities  are  considered  at  various 
locations  in  the  computational  model  of  the  mock  PC,  as  shown  in  Fig.  2-10.  The  various 
locations  considered  are  the  following: 

•  A,  in  the  power  supply  unit,  as  waves  entering  this  cavity  section  may  resonate  and 
may  generate  field  that  get  transmitted  to  every  other  location  in  the  PC  that  connects 
to  the  power  supply  unit. 

•  B,  directly  below  the  heat  sink  where  the  CPU  and  other  crucial  electronic  circuits 
may  be  situated  on  the  motherboard. 

•  C,  in  front  of  the  heat  sink  where  various  electronic  circuitry  is  located  on  the 
motherboard. 

•  D,  behind  the  CD  and  floppy  drive  where  the  ends  of  the  power  rails  from  the  power 
supply  unit  are  found. 
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Fig.  2-10.  Location  of  the  various  field  points  inside  the  moek  PC 


(a)  the  field  at  A 


(b)  the  field  at  B 

Fig.  2-11.  The  field  E^  at  (a)  A,  (b)  B,  (e)  C,  and  (d)  D  loeated  in  the  moek  PC  shown  in 

Fig.  2-10  (Continued) 
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(c)  the  field  at  C 


(d)  the  field  at  D 

Fig.  2-11.  The  field  E^  at  (a)  A,  (b)  B,  (e)  C,  and  (d)  D  loeated  in  the  moek  PC  shown  in 

Fig.  2-10 

Looking  at  the  field  E^  at  the  locations  mentioned  above,  one  can  summarize  the  results 
shown  in  Fig.  2-11. 

For  the  field  E^  at  point  A,  which  is  located  in  the  power  supply  unit,  the 

disturbance  appears  later  than  that  at  other  locations.  The  initial  amplitude  of  the  field  at 
this  location  is  fairly  small  compared  to  the  amplitude  of  the  fields  at  other  locations.  The 
low  initial  amplitude  indicates  very  little  transient  signal  penetration  into  the  power 
supply  unit  through  the  vent  slots.  As  time  passes,  the  fields  within  the  power  supply  unit 
start  bouncing  back  and  forth  between  the  walls  and  create  a  steady-state  maximum 
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amplitude  for  the  field,  whieh  ean  be  attributed  to  the  low  dissipation  of  fields 
through  the  vent  slots. 

For  the  field  at  point  B,  loeated  below  the  heat  sink,  there  is  good  penetration  of 
the  fields  through  the  slots  loeated  at  the  rear  of  the  PC  so  a  disturbanee  here  appears 
mueh  earlier  and  has  the  highest  initial  amplitude  of  all  observed  loeations.  In  faet,  the 
region  below  the  heat  sink  appears  to  fields  as  a  parallel  plate  waveguide,  so  the 
presenees  of  strong  fields  is  expeeted  as  is  the  lingering  of  the  fields  in  this  region  long 
after  the  souree  has  died.  The  existenee  of  sueh  high-amplitude  fields  where  the  CPU  and 
other  erueial  eleetronie  eireuits  may  be  situated  on  the  motherboard  may  result  in 
degradation  of  their  performanee. 

For  loeation  C,  in  front  of  the  heat  sink,  the  amplitude  of  the  field  overall  is  fairly 
small  as  the  heat  sink  aets  as  a  barrier  for  the  fields  penetrating  from  the  slots  loeated 
direetly  behind.  The  field  initially  has  small  amplitude  as  the  heat  sink  obstruets  the 
penetrating  transients  from  reaehing  the  loeation  of  the  field  point  and  then  inerease  to  a 
maximum  value,  which  can  be  attributed  to  the  good  penetration  of  the  fields  through  the 
slots  in  the  rear  of  the  PC.  The  field  then  reaches  a  minimum  after  the  source  dies  out. 
Finally,  the  field  increases  again  due  to  the  effect  of  fields  lingering  in  the  PC  long  after 
the  source  has  died. 

The  amplitude  of  the  field  at  D,  behind  the  CD  and  floppy  drive,  is  small 

compared  to  the  field  at  B.  The  field  has  high  initial  amplitude  owing  to  good 
penetration  of  the  transients  thought  the  vent  slots  and  decreases  gradually  to  the 
amplitude  of  the  fields  lingering  in  the  rectangular  chassis  of  the  mock  PC. 
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IV.  Effects  Of  The  Propagation  Path  On  Signae  Penetration 


To  further  analyze  the  effeets  of  the  propagation  path  on  the  penetration  of  signals 
into  an  eleetronie  system,  a  bent  wire  is  ineluded  in  the  moek  PC,  as  shown  in  Fig.  2-12, 
to  simulate  a  power  supply  eord  eoming  from  outside  into  the  system.  The  thin  wire  of 
radius  a  =  0.003m  emerges  from  the  ground  plane  at  a  distanee  5A  away  from  the  rear  of 
the  PC  and  rises  to  a  height  of  5 A  from  the  ground  plane.  The  wire  then  bends  and 
enters  the  PC  through  a  square  hole  3A  x  3A  in  size  and  then  splits  into  two  wires,  one  of 
whieh  goes  to  the  power  supply  unit  while  the  other  goes  to  a  dise  drive,  as  shown  in  Fig. 
2-12.  The  two  wires  after  the  split  in  the  bent  wire  are  made  to  run  very  elose  to  eaeh 
other,  as  shown  in  Fig.  2-12,  to  see  any  signifieant  effeets  due  to  propagation  path  and  the 
eoupling  effeet  between  elosely  spaeed  wire  bundles  in  a  real  PC.  The  numerieal 
proeedure  employed  in  modeling  the  bent-wire  strueture  is  deseribed  in  detail  in  Chapter 
4.  The  strueture  shown  in  Fig.  2-12  is  illuminated  with  a  thin-wire  monopole  of  height 
I  =  20A  and  radius  a  =  0.003m,  plaeed  15A  direetly  behind  the  PC  on  its  eenterline.  The 
monopole  is  fed  with  a  transmission  line  of  eharaeteristie  impedanee  Zg  =  50Q  and  is 

exeited  with  the  differentiated  Gaussian  pulse  shown  in  Fig.  2-3. 

To  view  the  effeets  of  the  propagation  path  (bent- wire  strueture)  on  the  transient 
signal  penetration  into  the  eleetronie  system,  eonsider  the  planes  along  the  y-dimension 
and  the  z-dimension  in  the  eross-seetion  ineluding  the  wire,  as  shown  in  Fig.  2-13.  With 
the  strueture  being  illuminated  by  the  monopole,  one  should  expeet  to  see  the  transient 
signal  penetration  inerease  as  the  fields  from  outside  get  eoupled  to  the  bent  wire 
strueture  leading  into  the  PC.  Indeed,  while  the  wire  strueture  enhanees  the  ability  of  a 
signal  to  reaeh  the  interior  of  the  PC  due  to  eonduetive  eoupling,  it  also  makes  it  possible 
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for  trapped  fields  lingering  inside  the  PC  to  escape  into  the  exterior  region,  which 


shows  up  at  later  time  steps  in  the  FDTD  simulation. 


Fig.  2-12.  The  computational  model  of  the  mock  PC  with  the  bent-wire  structure  to 

simulate  a  power  supply  cord 
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Fig.  2-13.  The  plane  along  the  y-direction  in  the  cross-section  of  the  wire  of  the  mock  PC 

with  the  bent-wire  structure 
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(a)  Time  step  =140 

Fig.  2-14.  The  field  in  the  plane  along  the  y-dimension  in  the  cross-section  of  the  wire 
shown  in  Fig.  2-13  for  (a)  time  step  =  140  and  (b)  time  step  =  1000  (Continued) 
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(b)  time  step  =  1000 

Fig.  2-14.  The  field  in  the  plane  along  the  y-dimension  in  the  cross-section  of  the  wire 
shown  in  Fig.  2-13  for  (a)  time  step  =  140  and  (b)  time  step  =  1000 


When  the  antenna  initially  illuminates  the  mock  PC  with  the  bent-wire  structure, 
fields  couple  to  the  bent  wire  and  make  their  way  into  the  PC,  as  shown  in  the  plot  in  Fig. 
2-14  (a).  In  Fig.  2-14  (b)  at  time  step  =  1000,  which  is  long  after  the  source  outside  has 
died  out,  one  can  still  see  the  fields  lingering  in  the  rectangular  chassis  of  the  PC,  but  the 
bent-wire  structure  clearly  aids  in  dissipating  the  fields  by  enhancing  radiation  out  the 
rear  of  the  PC. 

Now  consider  the  plane  along  the  z-dimension  in  the  cross-section  of  the  wire  as 
shown  in  Fig.  2-15. 
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Fig.  2-15.  The  plane  along  z-dimension  in  the  cross-seetion  of  the  wire  of  the  moek  PC 

with  the  bent  wire  strueture 


i  coordinate 

(a)  Time  step  =120 

Fig.  2-16.  The  field  in  the  plane  along  the  z-dimension  in  the  eross  section  of  the  wire 
shown  in  Fig.  2-15  for  (a)  time  step  =  120  and  (b)  time  step  =  600  (Continued) 
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(b)  Time  step  =600 


Fig.  2-16.  The  field  in  the  plane  along  the  z-dimension  in  the  cross  section  of  the  wire 


shown  in  Fig.  2-15  for  (a)  time  step  =  120  and  (b)  time  step  =  600 


The  plots  of  the  field  for  time  step  =  120  in  Fig.  2-16  (a)  illustrates  the  coupling  of 

the  spurious  electromagnetic  energy  to  the  wire  leading  into  the  PC.  From  Fig.  2-16  (b) 
one  can  see  the  fields  still  lingering  in  the  PC  long  after  the  source  outside  has  died  and 
also  the  dissipation  of  the  fields  brought  about  by  the  bent  wire.  Thus,  cables  running  into 
the  electronic  system  can  cause  electromagnetic  interference  that  affects  the  operation  of 
digital  circuits  embedded  deep  within  the  system.  Again,  to  obtain  a  better  understanding 
of  the  electromagnetic  waves  entering  an  electronic  system  from  an  outside  source,  one 
can  consider  the  field  intensities  at  various  locations  in  the  computational  model  of  the 
mock  PC  shown  in  Fig.  2-17. 
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Amplitude 


Fig.  2-17.  Location  of  the  various  field  points  inside  the  moek  PC  with  the  bent-wire 

structure 


(a)  the  field  at  A 


(b)  the  field  at  B 

Fig.  2-18.  The  field  E^  at  (a)  A,  (b)  B,  (c)  C,  and  (d)  D  loeated  in  the  moek  PC  with  the 
bent- wire  strueture  shown  in  Fig.  2-17  (Continued) 
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(c)  the  field  at  C 


(d)  the  field  at  D 

Fig.  2-18.  The  field  E^  at  (a)  A,  (b)  B,  (c)  C,  and  (d)  D  located  in  the  mock  PC  with  the 

bent- wire  structure  shown  in  Fig.  2-17 

Comparing  the  field  E^  at  various  locations  in  Fig.  2-10  for  the  cases  without  the 

bent-wire  structure,  as  in  Fig.  2-11,  and  with  the  bent-wire  structure,  as  in  Fig.  2-18,  one 
can  summarize  the  results  in  Fig.  2-18  as  follows. 

For  point  A  in  the  power  supply  unit  the  field  at  this  location  enters  at 
approximately  the  same  time  as  the  fields  in  other  locations  as  the  transient  fields  now  get 
coupled  to  the  wire  and  enter  the  PC.  The  initial  amplitude  of  the  field  is  fairly  small  but 
higher  than  the  results  of  Fig.  2-11  (a).  Again,  as  time  passes  the  penetrated  fields  start 
bouncing  back  and  forth  between  the  walls  of  the  power  supply  unit.  This  is  observed 
much  earlier  than  in  the  results  for  the  model  without  the  bent  wire,  which  are  shown  in 
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Fig.  2-11  (a).  The  amplitude  of  the  field  at  this  loeation  gradually  deereases  for  late 
time,  and  this  ean  be  attributed  to  earrying  off  of  energy  from  the  unit  via  the  bent-wire 
strueture. 

Point  B  is  below  the  heat  sink  so  the  field  eoupling  to  the  bent-wire  strueture 
leading  results  in  a  field  with  higher  amplitude  than  that  of  field  in  Fig.  2-11  (b).  Due 

to  the  lingering  of  the  fields  inside  the  ehassis  of  the  PC  long  after  the  souree  has  died  out 
and  the  wire  strueture  field  eoupling  effeet,  the  field  at  this  loeation  remains  fairly  high 
even  for  late  time. 

In  front  of  the  heat  sink,  at  point  C,  the  penetrated  transient  signal  is  very  similar 
to  Fig.  2-11  (e)  with  the  slightly  higher  amplitude  attributed  to  enhaneed  eoupling  due  to 
the  bent-wire  strueture. 

Behind  the  CD  and  floppy  drive  at  D  the  field  has  very  high  amplitude  when 
eompared  to  Fig.  2-11  (d),  thus  showing  the  effeets  of  propagation  path  on  the 
penetration  of  signals  into  the  eleetronie  system.  The  field  remains  high  for  long  periods 
of  after  the  souree  has  died  out  due  to  the  eoupling  of  fields  between  the  elosely  spaeed 
wires  leading  to  the  power  supply  unit  and  the  CD  drive  and  the  field  seattering  around 
the  sharp  bends. 


V.  Summary 

This  ehapter  presents  an  overview  of  eleetromagnetie  interferenee  effeets  relating 
to  the  performanee  of  an  eleetronie  system.  The  effeets  of  transient  eleetromagnetie 
waves  on  highly  sensitive  digital  eireuits  in  eleetronie  systems  are  very  eomplex  and 
diffieult  to  understand  as  simple  mismatehes,  slots,  or  diseontinuities  ean  eause  eoupling 
in  the  inlaid  eireuits.  The  results  presented  for  an  eleotromagnetieally  replieated  PC 
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model  helps  one  to  prediet  the  properties  of  transient  signals  as  they  make  their  way 
through  a  propagation  path  from  the  exterior  of  a  strueture  to  the  terminals  of  a  deeply 
embedded  digital  system  or  circuit  and  to  catalog  the  properties  of  signals  that  might  be 
expected  to  reach  various  locations  in  the  electronic  system.  Observing  the  dominant 
field  component  at  certain  specific  locations  in  the  mock  PC  gives  one  a  better 
understanding  of  the  behavioral  patterns  of  the  penetrating  fields  at  these  locations.  Metal 
wires  both  with  and  without  bends  are  likely  to  have  reactive  impedances  that  degrade  the 
performance  of  complex,  high-speed  circuitry  surrounding  them.  The  model  of  the  mock 
PC  with  the  bent-wire  structure  showed  the  effects  of  the  wires  in  the  electronic  systems. 
The  most  significant  feature  observed  in  the  case  of  adding  the  bent-wire  structure  is  the 
enhancement  of  transient  signal  penetration  into  the  system,  attributed  to  the  fields 
getting  coupled  to  the  wire  from  outside.  These  effects  of  transient  signal  penetration  into 
electronic  systems  can  be  further  studied  by  adding  more  realistic  detail  to  determine  the 
causes  and  effects  of  specific  components  of  interest  and  by  considering  other  structures 
and  excitations  (pulses,  pulse  trains,  etc.). 
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Chapter  3 

Extrapolation  In  Time  And  Frequency  Domains 


Using  Orthonormal  Associate  Hermite 

Functions 

I.  Introduction 

The  effects  of  the  transmission  path  and  environment  on  the  characteristics  of  a 
transient  signal  due  to  a  high-power  microwave  (HPM)  source  as  it  passes  through  an 
electronic  system  is  of  importance  in  discovering  how  the  operation  of  digital  circuits 
might  be  altered,  since  such  signals  are  unexpected  and  contain  spurious  electromagnetic 
energy.  If  a  signal  of  a  specified  form  enters  a  system  at  a  given  point  it  is  useful  to  know 
the  salient  features  of  the  resulting  signal  that  reaches  a  location  within  the  system  where 
a  susceptible  digital  circuit  is  located.  Moreover,  it  is  useful  to  know  which  features  of 
the  induced  signal  are  primarily  those  peculiar  to  the  entering  signal  and  which  features 
are  more  influenced  by  the  properties  of  the  transmission  environment.  Addressing  these 
issues  should  help  one  gain  an  appreciation  for  the  nature  of  a  spurious  signal  arriving  at 
the  input  of  a  digital  circuit  embedded  deep  within  a  complex  system,  and  perhaps  that 
information  can  be  used  to  eliminate  deleterious  effects  on  the  operation  of  digital 
circuits. 

With  the  powerful  Maxwell  equations  solvers  available  today  to  perform 
electromagnetic  simulation,  it  is  reasonable  to  assume  that  they  can  be  useful  to 
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characterize  sueh  effeets.  However,  given  the  eomplexity  of  real-world  systems  this 
still  remains  a  daunting  task,  espeeially  when  wideband  charaeterization  is  needed.  For 
example,  the  use  of  the  popular  FDTD  method  for  modeling  eleetronics  inside  enclosures 
imposes  serious  memory,  time,  and  aceuraey  constraints  in  generating  the  late-time 
response.  Also,  the  use  of  MOM  or  FEM  when  generating  the  high-frequency  response 
ean  be  eomputationally  overwhelming.  Reeently,  researehers  have  developed  wideband 
extrapolation  techniques  that  are  often  capable  of  overeoming  some  of  these  limitations 
by  using  only  early-time  and  low-frequeney  data  [1].  By  using  only  the  early  time  and 
low  frequency  data,  which  contain  mutually  complementary  information,  one  ean 
simultaneously  extrapolate  the  entire  system  response  in  both  domains  [2]. 

The  focus  of  this  chapter  is  to  examine  how  well  these  extrapolation  techniques 
work.  The  FDTD  method  is  used  to  obtain  the  early-time  response,  and  a  MOM-based 
(any  frequeney  domain  MOM  or  FEM  code  could  be  used)  is  used  to  obtain  the  low- 
frequeney  response.  Attempts  are  made  to  aseertain  the  aeeuraey  of  the  approaeh  as  well 
as  the  CPEl  time  savings  measured  against  using  either  the  EDTD  or  MOM  tools 
separately. 


11.  Orthonormal  Associate  Hermite  Functions 

Using  only  the  early  time  and  low  frequency  data,  which  contain  mutually 
complementary  information,  one  can  simultaneously  extrapolate  the  entire  system 
response  in  both  domains  using  the  orthonormal  associate  hermite  (AH)  functions.  The 
isomorphism  between  an  AH  function  and  its  Eourier  transform  allows  one  to  work  in 
both  time  and  frequency  domains  simultaneously  [1-3]. 
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The  associate  hermite  functions  (x)  are  defined  through  hermite 


polynomials  (x)  as 


2 


where  n  is  order  of  the  hermite  polynomials.  The  hermite  polynomials  are  recursively 
expressed  as 


«oW  =  l 

//l(x)  =  2x  (8.2) 

Hn  (^)  =  (x)  -2(n  - \)H^_2  (x),  for  n  >  2 


using  (8.2)  the  AH  functions  can  be  expressed  as 

2 

(8.3) 

for  n  >  2 

The  plot  of  AH  functions  for  orders  zero  to  four  are  as  shown  in  Fig.  3-1.  The  AH 
functions  are  also  known  as  orthonormal  AH  functions  because  their  polynomials  are 
mutually  orthogonal  and  form  a  complete  set  in  the  interval  [-go,-i-go]  [2],  given  as 
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(8.4) 


I  hi  (x) hj  (x)  Jx  =  dy 


where  Stj  is  the  Kronecker  delta  for  which 


dy  =  1,  for  i  =  j 

5y  =  0,  otherwise 


(8.5) 


Fig.  3-1  Associate  hermite  functions  of  orders  zero  to  four 


As  shown  in  Fig.  3-1,  the  AFl  polynomials  of  even  order  are  even  functions  while  those 
of  odd  order  are  odd  functions 


^2m  (^)  =  hm  (-x) 

^2m+l(^)  =  -^2m+l(-^) 

for  any  value  of  m.  Finally,  if  / (x)  is  a  finite  interval  piecewise  smooth  function  in  the 
interval  [-k,A:]  and 
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(8.7) 


/  [x)dx<cc 


then  f{x)  can  be  represented  as 


f{x)=^a„h„  (x),  for  —  oo<x<oo 

n=0 


with 


(8.8) 


a„  =\^^f{x)h„{x)dx 


(8.9) 


This  series  exhibits  pointwise  convergence  to  f{x)  wherever  f{x)  is  continuous  and  to 


/U  +/U 


12  at  points  of  discontinuity.  From  the  property  of  isomorphism,  the 

AH  functions  are  eigenfunctions  of  the  Fourier  transform  operator.  To  clarity,  if 
is  defined  to  be  the  Fourier  transform  of  (x)  then 


(A)  =  ^_^h„  {x)Qxp{-j27rXx)dx  =  {-jf  h„  (X)  (8.10) 

where  j  =  .  Hence,  if  F(X )  is  considered  to  be  the  Fourier  transform  of  / (x) ,  then 

one  can  write 

F{X)  =  \"j{xyx^{-j2nXx)  (8.11) 

Substituting  (8.8)  into  (8.11) 
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(8.12) 


and  hence  from  (8.10) 


F{X)=J,a„{-j)"  h„{X)  (8.13) 

n 

From  (8.8)  and  (8.13),  one  can  conclude  that  the  same  hermite  functional  representation 
can  be  used  to  approximate  both  time-domain  and  frequency-domain  electromagnetic 
responses. 


III.  The  Mathematical  Formulation 

Let  /(t)  be  the  time  domain  electromagnetic  response  due  to  an  incident  stimulus 
(t) .  Using  the  orthonormal  AH  functions,  one  can  represent  the  time-domain 
electromagnetic  response  as 


N 


n=0 


(8.14) 


where  is  a  scaling  factor  for  the  time  variable.  As  in  (8.8),  the  time-domain  response 
obtained  using  the  time-domain  techniques  cannot  go  to  infinity  but  only  up  to  some 
value  Nf  because  one  is  only  interested  in  the  obtaining  the  early-time  response. 
Similarly,  the  frequency-domain  response  can  also  be  represented  in  terms  of  AH 
functions.  From  (8.13)  the  frequency-domain  response  F’(/)  can  be  written  as 
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(8.15) 


Again,  as  the  frequency-domain  response  computed  is  finite  with  N j-  points.  Equation 
(8.15)  can  further  be  expressed  as 


F{f)=^{-jT^K{f/s2) 

n=0  V^2 

^hnif  !  Si)-  j^^f^hn+lif  !  Si) 
^Jsi  ^jsi 

F{f)  =  Fi,{f)  +  jFj{f) 


NH 


n=0 


(8.16) 

(8.17) 

(8.18) 


Again  $2  is  the  scale  factor  for  frequency.  The  scale  factors  in  either  domain  are  related 
as 


1 

^  [Itisi) 


(8.19) 


This  extrapolation  technique  is  most  effective  when  the  causal  time  domain  signal  f{t), 
to  be  expanded  in  terms  of  AH  functions,  is  centered  about  t  =  t^,  where  is  around  half 
the  time  support  of  /  (t)  [1-2]  because  the  AH  functions  provide  equal  support  on  either 

side.  Centering  the  expansion  would  require  fewer  terms,  hence  the  corresponding 
transform  in  the  frequency  domain  is  given  as 

f{t  +  tQ)  =  F{f)QXTp{j27iftQ)  (8.20) 
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The  scaling  factors  are  crucial  because  they  determine  the  amount  of  support 
provided  by  the  AH  functions  to  the  time  and  frequency  domain.  Using  sufficient 
amounts  of  low-frequency  and  early  time-domain  data,  a  properly  chosen  order  of 
expansion  N,  and  well  chosen  scale  factors  and  S2 ,  one  can  obtain  the  missing  data 

(i.e.,  the  high-frequency  and  late  time-domain  data).  The  order  of  expansion  N  is  obtained 
by  trial  and  error,  by  using  a  large  value  initially  and  choosing  a  cutoff  for  the  magnitude 
of  the  coefficients  and  discarding  the  ones  that  die  out.  An  important  point  to  be 
considered  while  choosing  the  order  is  that  an  unnecessarily  large  N  introduces 
oscillations  in  the  results. 


IV.  Matrix  Formulation 

The  matrix  representation  of  the  early  time-domain  data  in  terms  of  the  hermite 
functions  (8.14)  is  expressed  as 


/2o(hAi)  ^i(hAi) 
K{h/s\)  hi{t2lsi)  ■■ 

% 

=  ^/N 

/(h) 

^7V-l(^7V, /■^i) 

N,xN 

Nxl 

/('«,) 

The  real  part  of  the  low  frequency  data  A(/)  from  (8.16)  can  be  represented  in  matrix 
form  as 
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-hifilsi)  •••  (-1)' 

-hifihi)  •••  (“O'  ^ 


-h 


2\JN 


4^ 


FrUi) 

pRifl) 


(8.22) 


Ft 


R\JN 


J7V,xl 


Similarly,  the  imaginary  part  of  the  low  frequency  data  f^(/)  from  (8.18)  is 


-hiifihi) 

hifihi)  •• 

(“O' %-2 (/1A2) 

hifihi)  ■■ 

~h  {^fNf 

1 

1 _ 

=  44 

^  F, (A)  ^ 
Flifl) 

fxl 

1 

_ _ _ 1 

(8.23) 


Combining  (8.21),  (8.22),  and  (8.23), 
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/lo(hM) 

hihhi) 

^N-iihhi) 

^0  (%,  /■^l ) 

/■^i) 

-hifihi) 

(“O'  ^^7V-i(/iA2) 

k)  [In^ 

~^2  [/n^  1^2) 

(“O'  ^  ^7V-1  (/tv,. /■^2) 

hifihi) 

(“O'  %-2  (/l  A2) 

~h  (/tv,  j^i) 

h  (/tv^.  1^2) 

(“O'  %-2  (/tv, /■^2) 

l(Af,+27VjxAf 


yfllfih) 


aj 


(/l) 


«2 

a]^_l 


Nx\ 


^S'^■;(/l) 


-l(7V,+27V^)xl 


(8.24) 


Therefore,  in  (8.24)  using  the  early  time-domain  data  from  /(h), /(t2 ),•••,/ )  and 


the  low  frequency-domain  data  F  (/),  F  (/2  ),...,  F  ( /^  |  the  unknown  coefficients  of 
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the  expansion  are  obtained  by  solving  a  least-squares  problem  using  singular  value 
deeomposition  [1-3].  Using  the  obtained  unknown  eoeffieients  one  ean  then  eonstruet 
funetion  representations  to  extrapolate  both  early  time  and  low  frequency  domain  data  to 
late-time  and  high-frequency  domain  data. 


V.  Results 

To  validate  the  technique  of  using  only  early-time  and  low-frequency  data,  which 
contain  mutually  complementary  information,  to  simultaneously  extrapolate  the  entire 
system  response  in  both  domains  using  the  orthonormal  associate  hermite  (AH) 
functions,  consider  a  simple  dipole  antenna  of  length  /  =  1.0m  and  radius  a  =  .001m 
loaded  with  a  resistive  load  of  i?  =  lOOQ  in  free  space.  The  current  in  the  load  of  this 
dipole  is  characterized  in  the  time  domain  using  FDTD  and  in  the  frequency  domain 
using  MOM,  as  described  in  Chapter  5.  The  number  of  unknowns  used  to  approximate 
the  current  distribution  on  the  antenna  for  MOM  is  63,  and  the  FDTD  cell  size  used  to 
obtain  the  time-domain  results  is  Ax  =  Ay  =  Az  =  A  =  l/63m.  The  load  is  positioned  at 
the  48**^  segment,  which  corresponds  to  the  middle  of  the  monopole.  One  can  use  any 
arbitrary  excitation  as  a  voltage  source  at  the  feed  to  generate  the  current  in  the  load,  but 
the  excitation  used  here  is  a  Gaussian  source  expressed  as 


(8.25) 


where  /  =  [t-tQ) j a ,  tg  is  a  delay  used  to  make  the  pulse  rise  smoothly  from  zero  for 
time  t  <  0  to  its  value  at  t,  and  cr  represents  the  standard  deviation  of  the  Gaussian 
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distribution.  The  Gaussian  excitation  waveform  and  its  frequency  spectrum  are 


shown  in  Fig.  3-2. 


(c)  Gaussian  excitation  waveform  in  the  time  domain 


(d)  The  frequency  spectrum  of  the  excitation  waveform  of  (a) 

Fig.  3-2.  The  Gaussian  voltage  source  introduced  at  the  feed  point  of  the  thin-wire 

antenna 


To  obtain  the  frequency  domain  response  of  the  loaded  dipole  subjected  to  the  above 
Gaussian  pulse,  the  frequency  domain  data  obtained  using  the  MOM  has  to  be 
multiplied  by  the  spectrum  of  the  time-domain  Gaussian  pulse. 

The  time-domain  results  shown  in  Fig.  3-3  are  obtained  from  ^  =  0  to 
^  =  13.4502 ns  (500  data  points)  using  FDTD,  and  the  frequency-domain  results 
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shown  in  Fig.  3-4  (a)  are  obtained  from  DC  to  /  =  3. 9980 GHz  (2,000  data  points) 
using  MOM.  The  frequency  response  of  the  wire  antenna  in  Fig.  3-4  (b)  is  computed 
by  multiplying  the  frequency-domain  data  obtained  by  MOM  with  the  spectrum  of 
the  Gaussian  pulse  of  Fig.  3-2  (b). 


Fig.  3-3.  The  time-domain  data  for  the  current  in  the  load  of  a  dipole  antenna  of  length 
/  =  1.0  m  and  radius  a  =  .001m  loaded  with  a  resistive  load  of  i?  =  lOOQ  in  free  space 


(a)  Frequency  domain  data  obtained  via  MOM 


Fig.  3-4.  The  frequency-domain  data  for  the  current  in  the  load  of  a  dipole  antenna  of 
length  /  =  1.0  m  and  radius  a  =  .001m  loaded  with  a  resistive  load  of  i?  =  lOOQ  in 

free  space  (Continued) 
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(b)  Frequency  response  of  the  system 


Fig.  3-4.  The  frequency-domain  data  for  the  current  in  the  load  of  a  dipole  antenna  of 
length  I  =  1.0  m  and  radius  a  =  .001m  loaded  with  a  resistive  load  of  i?  =  lOOQ  in 

free  space 

As  mentioned  earlier,  one  can  extrapolate  the  time-  and  frequency-domain  data  by 
using  sufficient  low-frequency  and  early  time-domain  data  {Nf),  a  properly 

chosen  order  of  expansion  N,  and  scale  factors  and  S2-  For  the  time  being,  consider  an 
expansion  of  order  A  =  100  having  scale  factors  chosen  so  that  the  frequency  range  of  the 
available  data  is  mapped  to  (-6,  6) .  Assuming  that  only  the  first  180  data  points  are 

available  in  the  time  domain  (A^)  and  only  the  first  700  data  points  are  available  in  the 
frequency  domain  the  matrix  equation  given  by  (8.24),  when  solved  for  the 

unknown  coefficients,  can  be  used  to  extrapolate  the  time-domain  data  to  t  =  13.4502ns 
(500  data  points)  and  the  frequency-domain  data  to  /  =  3.9980 GHz  (2,000  data  points), 
as  shown  in  Fig.  3-5. 
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(a)  Time-domain  response 


Frequency  steps 


(c)  Real  part  of  the  frequency  response 


Frequency  steps 


(d)  Imaginary  part  of  the  frequency  response 


Fig.  3-5.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequency-domain  response,  (c)  imaginary  part  of  frequency-domain  response,  and  (d) 
the  unknown  expansion  coefficients  of  the  AFIF  for  a  resistively  loaded  dipole  of 
length  I  =  1.0m,  a  =  .001m,  and  R  =  lOOQ  for  Nf  =  180,  N =  700  ,  N  =100, 

and  the  frequency  range  of  the  available  data  mapped  to  (-6,  6)  (Continued) 
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(e)  The  unknown  coefficients 


Fig.  3-5.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequency-domain  response,  (c)  imaginary  part  of  frequency-domain  response,  and  (d) 
the  unknown  expansion  coefficients  of  the  AFIF  for  a  resistively  loaded  dipole  of 
length  /  =  1.0m,  a  =  .001m,  and  i?  =  100Q  for  Nj  =180,  N  =  700,  N  =100, 

and  the  frequency  range  of  the  available  data  mapped  to  (-6,  6) 


The  time-domain  waveform  is  centered  about  =2.5555(u5),  which  is  denoted 

by  an  asterisk  in  Fig.  3-5  (a).  The  data  assumed  to  be  available  in  both  domains  is 
separated  from  the  extrapolated  data  by  vertical  dashed  lines.  The  results  are  incorrect 
and  show  significant/large  late  time  and  high  frequency  oscillations,  indicating 
insufficient  amount  of  support  being  provided  by  the  AH  functions  to  the  time  and 
frequency  domain  data.  This  can  be  corrected  by  increasing  the  scale  factor  so  that  the 
frequency  range  of  the  available  data  is  mapped  to  (-13,  13).  Given  the  same  time- 

domain  and  frequency-domain  data  as  before,  the  solution  of  the  matrix  equation  given 
by  (8.24)  for  the  unknown  coefficients  yields  the  results  shown  in  Fig.  3-6. 
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(a)  Time-domain  response 


(b)  Real  part  of  the  frequeney  response 


(e)  Imaginary  part  of  the  frequeney  response 


Fig.  3-6.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequeney-domain  response,  (e)  imaginary  part  of  frequency-domain  response  and  (d) 
the  unknown  expansion  coefficients  of  the  AHF  for  a  resistively  loaded  dipole  of 
length  /  =  1.0m,  a  =  .001m,  and  i?  =  100Q  for  Nf  =180,  N  f  =  700,  N  =100, 

and  the  frequency  range  of  the  available  data  mapped  to  (-13,  13)  (Continued) 
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(d)  The  unknown  coefficients 


Fig.  3-6.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequency-domain  response,  (c)  imaginary  part  of  frequency-domain  response  and  (d) 
the  unknown  expansion  coefficients  of  the  AFIF  for  a  resistively  loaded  dipole  of 
length  /  =  1.0m,  a  =  .001m,  and  i?  =  100Q  for  Nj  =180,  N  =  700,  N  =100, 

and  the  frequency  range  of  the  available  data  mapped  to  (-13, 13) 


The  time-domain  waveform  is  again  centered  about  =  2.5555  (n^)  and  denoted 

by  an  asterisk  in  Fig.  3-6  (a).  These  results  are  much  better  than  the  results  in  Fig.  3-5  but 
still  do  not  agree  with  the  original  data  even  in  the  range  of  the  available  data.  The 
amplitude  of  the  plot  of  unknown  coefficients  in  Fig.  3-6  (d)  tends  to  increase  for  higher 
orders  of  expansion,  indicating  that  an  insufficient  amount  of  either  time-  or  frequency- 
domain  data  is  considered.  Hence,  increasing  the  number  of  available  data  points  in  the 

frequency  domain  (^  Ay  ^  from  700  to  1000  data  points  and  using  the  same  order  of 

expansion,  scale  factor,  and  number  of  time-domain  data  points  as  before,  much  more 
agreeable  results  are  obtained,  as  shown  in  Fig.  3-7. 
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—  Original  data 


(a)  Time-domain  response 


(b)  Real  part  of  the  frequency  response 


(c)  Imaginary  part  of  the  frequency  response 


Fig.  3-7.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequency-domain  response,  (c)  imaginary  part  of  frequency-domain  response,  and  (d) 
the  unknown  expansion  coefficients  of  the  AHF  for  a  resistively  loaded  dipole  of 
length  /  =  1.0m ,  a  =  .001m,  and  R  =  lOOQ  for  Nf  =  180,  N f  =  1000 ,  N  =100, 

and  the  frequency  range  of  the  available  data  mapped  to  (-13,  13)  (Continued) 
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Number  of  unknowns 


(d)  The  unknown  coefficients 


Fig.  3-7.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequency-domain  response,  (c)  imaginary  part  of  frequency-domain  response,  and  (d) 
the  unknown  expansion  coefficients  of  the  AFIF  for  a  resistively  loaded  dipole  of 
length  /  =  1.0m ,  a  =  .001m ,  and  R  =  lOOQ  for  =  180 ,  =  1000 ,  N  =100, 

and  the  frequency  range  of  the  available  data  mapped  to  (-13, 13) 


The  time-domain  waveform  is  centered  about  =  2.5555  (n^)  ,  denoted  by  an  asterisk  in 

Fig.  3-7  (a),  and  the  data  assumed  to  be  available  in  both  domains  is  separated  from  the 
extrapolated  data  by  vertical  dashed  lines.  The  plot  of  the  unknown  coefficients  shows 
the  coefficients  dying  out  for  higher  orders  of  expansion,  which  is  expected  [1-3].  From 
Fig.  3-7  one  can  say  that  the  extrapolated  data  is  agreeable  in  comparison  to  the  original 
data.  The  slight  mismatch  between  the  original  data  and  extrapolated  data  in  Fig.  3-5  (a) 
for  later  time  steps  can  be  attributed  to  grid  dispersion  errors  in  FDTD,  which  are 
explained  in  Chapters  2  and  5.  Hence,  for  a  given  time -bandwidth  product  of  9.6207,  one 
is  able  to  extrapolate  to  a  time -bandwidth  product  of  53.6662. 

The  CPU  times  taken  to  compute  the  assumed  available  time-domain  and 
frequency-domain  data  on  a  2.08  GHz  processor  with  2.0  GB  of  RAM  are  13.47  and 
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23.07  minutes,  respectively,  yielding  a  total  time  of  36.54  minutes.  If  one  were  to 
compute  the  entire  range  of  time-domain  and  frequency-domain  data  without  using 
extrapolation  techniques,  it  would  take  37.30  minutes  for  time-domain  calculations  and 
45.26  minutes  for  frequency-domain  calculations.  The  amount  of  CPU  time  saved  for  the 
simple  case  of  a  loaded  dipole  antenna  is  46.02  minutes.  This  CPU  time  savings  becomes 
even  more  significant  when  calculating  the  electromagnetic  responses  of  complex  cavity 
structures  in  both  domains,  a  process  that  is  normally  very  computationally  intensive  and 
time-consuming.  It  is  clear  that  the  extrapolation  techniques  can  significantly  reduce  the 
amount  of  CPU  time  used. 

Now  consider  a  simple  dipole  antenna  of  length  /  =  2.0m  and  radius  a  =  .005m 
in  free  space.  The  current  at  the  feed  point  of  this  dipole  is  characterized  in  the  time 
domain  using  FDTD  and  in  the  frequency  domain  using  MOM,  as  described  in  Chapter  5. 
The  number  of  unknowns  used  to  approximate  the  current  distribution  on  the  antenna  for 
MOM  is  601,  and  the  FDTD  cell  size  used  to  obtain  the  time-domain  results  is 
Ax  =  Ay  =  Az  =  A  =  0.0392  m .  The  time-domain  results  are  obtained  from  t  =  0  to 
t  =  19.6214ns  (300  data  points)  using  FDTD,  and  the  frequency-domain  results  shown  in 
Fig.  3-8  (a)  are  obtained  from  DC  to  /  =  4.0GHz  (2,000  data  points)  using  MOM.  The 
frequency  response  of  the  wire  antenna  is  again  obtained  by  multiplying  the  frequency- 
domain  data  obtained  by  MOM  with  the  spectrum  of  the  incident  voltage  waveform  used 
in  the  FDTD  code.  Assuming  that  only  the  first  100  data  points  are  available  in  the  time 
domain  and  only  the  first  650  data  points  are  available  in  frequency  domain,  one  can 
solve  the  matrix  equation  given  by  (8.24)  for  the  unknown  coefficients  and  successfully 
extrapolate  the  time-domain  data  to  i  =  19.6214  ns  (300  data  points)  and  the  frequency- 
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domain  data  to  /  =  4.0GHz  (2,000  data  points),  as  shown  in  Fig.  3-8.  The  order  of 


the  expansion  is  chosen  to  be  =  90  and  the  plot  of  the  unknown  coefficients  of  (9.24) 
is  as  shown  in  Fig.  3-8  (d).  The  time-domain  waveform  is  centered  about 
=4.9054(^5),  which  is  denoted  by  an  asterisk  in  Fig.  3-8  (a).  The  choice  of  I2  was 

made  so  that  the  frequency  range  of  the  available  data  is  mapped  to  (-12,  12)  . 


(a)  Time-domain  response 


(b)  Real  part  of  the  frequency  response 


Fig.  3-8.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequency-domain  response,  (c)  imaginary  part  of  frequency-domain  response,  and 
(d)  the  unknown  expansion  coefficients  of  the  AHF  for  the  current  at  the  feed  point 
of  a  dipole  /  =  2.0  m  and  a  =  .005  m  for  =  100,  iVy-  =  650 ,  N  =90,  and  the 

frequency  range  of  the  available  data  mapped  to  (-12, 12)  (Continued) 
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(c)  Imaginary  part  of  the  frequency  response 


(d)  The  unknown  coefficients 


Fig.  3-8.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequency-domain  response,  (c)  imaginary  part  of  frequency-domain  response,  and 
(d)  the  unknown  expansion  coefficients  of  the  AHF  for  the  current  at  the  feed  point 
of  a  dipole  /  =  2.0  m  and  a  =  .005  m  for  =  100 ,  N  =  650  ,  N  =90,  and  the 

frequency  range  of  the  available  data  mapped  to  (-12, 12) 


In  Fig.  3-8  the  data  assumed  to  be  available  in  both  domains  is  separated  from  the 
extrapolated  data  by  vertical  dashed  lines.  The  extrapolated  data  is  agreeable  in 
comparison  to  the  original  data  with  a  slight  mismatch  between  the  original  data  and 
extrapolated  data  for  later  time  steps.  Hence,  for  a  given  time-bandwidth  product  of 
8.5026,  one  is  able  to  extrapolate  to  a  time-bandwidth  product  of  78.4857. 
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Lastly,  consider  another  case  of  a  simple  dipole  antenna  of  length  /  =  1.0m 
and  radius  a  =  .001m  in  free  spaee.  The  eurrent  at  the  feed  point  of  this  dipole  is  again 
eharaeterized  in  the  time  domain  using  FDTD  and  in  the  frequency  domain  using  MOM, 
as  described  in  Chapter  5.  The  number  of  unknowns  used  to  approximate  the  eurrent 
distribution  on  the  antenna  for  MOM  is  63,  and  the  FDTD  eell  size  used  to  obtain  the 
time-domain  results  is  Ax  =  Ay  =  Az  =  A  =  1/63  m  .  The  time-domain  results  are  obtained 
from  t  =  0  to  t  =  19.8550ns  (750  data  points)  using  FDTD,  and  the  frequeney-domain 
results  shown  in  Fig.  3-9  (a)  are  obtained  from  DC  to  /  =  4.0GHz  (2,000  data  points) 
using  MOM.  The  frequeney  response  of  the  wire  antenna  is  again  obtained  by 
multiplying  the  frequeney-domain  data  obtained  by  MOM  with  the  speetrum  of  the 
incident  voltage  waveform  used  in  the  FDTD  eode.  Assuming  that  only  the  first  190  data 
points  are  available  in  the  time  domain  and  only  the  first  900  data  points  are  available  in 
frequency  domain,  one  ean  solve  the  matrix  equation  given  by  (8.24)  for  the  unknown 
eoeffieients  and  suceessfully  extrapolate  the  time-domain  data  to  t  =  19.8550ns  (750 
data  points)  and  the  frequency-domain  data  to  /  =  4.0  GHz  (2,000  data  points),  as  shown 
in  Fig.  3-9.  The  order  of  the  expansion  is  chosen  to  be  A  =  100  and  the  plot  of  the 
unknown  eoeffieients  of  (9.24)  is  as  shown  in  Fig.  3-8  (d).  The  time-domain  waveform  is 
centered  about  C  =2.6473(n5),  which  is  denoted  by  an  asterisk  in  Fig.  3-9  (a).  The 

choiee  of  I2  was  made  so  that  the  frequeney  range  of  the  available  data  is  mapped  to 
(-13,13). 
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(a)  Time-domain  response 


(b)  Real  part  of  the  frequency  response 


(c)  Imaginary  part  of  the  frequency  response 


Fig.  3-9.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequency-domain  response,  (c)  imaginary  part  of  frequency-domain  response,  and 
(d)  the  unknown  expansion  coefficients  of  the  AHF  for  the  current  at  the  feed  point 
of  a  dipole  /  =  1.0m  and  a  =  .001m  for  =  190 ,  N  f  =  900  ,  N  =100,  and  the 

frequency  range  of  the  available  data  mapped  to  (-13,  13)  (Continued) 
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(d)  The  unknown  coefficients 


Fig.  3-9.  The  original  data  and  the  extrapolated  data  for  (a)  time-domain,  (b)  real  part  of 
frequency-domain  response,  (c)  imaginary  part  of  frequency-domain  response,  and 
(d)  the  unknown  expansion  coefficients  of  the  AHF  for  the  current  at  the  feed  point 
of  a  dipole  /  =  1.0m  and  a  =  .001m  for  =  190,  N =  900  ,  N  =100,  and  the 

frequency  range  of  the  available  data  mapped  to  (-13,  13) 


The  data  assumed  to  be  available  in  both  domains  is  separated  from  the 
extrapolated  data  by  vertical  dashed  lines  in  Fig.  3-9.  The  extrapolated  data  is  agreeable 
in  comparison  to  the  original  data  with  a  slight  mismatch  between  the  original  data  and 
extrapolated  data  for  later  time  steps.  Hence,  for  a  given  time-bandwidth  product  of 
9.0539,  one  is  able  to  extrapolate  to  a  time-bandwidth  product  of  79.4200. 

VI.  Summary 

In  this  chapter,  simultaneous  extrapolation  technique  for  both  time-domain  and 
frequency-domain  data  using  associate  hermite  functions  is  elucidated.  The  theoretical 
formulation  is  explained  in  detail,  and  results  are  clearly  illustrated.  The  results  for  a 
loaded  dipole  show  good  agreement  between  the  original  data  and  the  extrapolated  data, 
proving  that  these  techniques  work.  This  technique  requires  significantly  less  time  and 
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computer  memory  than  the  conventional  approaeh  involving  use  of  either  FDTD  or 


MOM  tools  separately. 
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